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We calculate the electromagnetic contribution to the pion mass difference, Am^ = m 2 K+ — m 2 , 
in the chiral limit through the VV — AA type vacuum polarization using Das-Guralnik-Mathur- 
Low-Young (DGMLY) sum rule. The calculation is made with two-flavors of dynamical over- 
lap fermions on a 16 3 x 32 lattice at a ^0.12 fm. The exact chiral symmetry of the over- 
lap fermion is essential to control the systematic error in the difference VV — AA. We obtain 
Anij = 1024(100) Me V 2 combining the lattice data with the perturbative contribution in the high 
momentum region evaluated by the operator product expansion. By analyzing the momentum de- 
pendence of the vacuum polarization, we also obtain pion decay constant f n and the low-energy 
constants L r l0 in the chiral limit. 
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1. Introduction 

The mass difference between charged pion (tt + ) and neutral pion (71°) is considered to be dom- 
inated by the electromagnetic (EM) contribution, which contains non-perturbative physics through 
the (off-shell) pion-pion-photon vertex. Determination of this quantity with good precision is nec- 
essary to determine up and down quark masses, which can be done in principle using lattice QCD. 

In 1967, Das, Guralnik, Mathur, Low and Young [[lj] derived a sum rule, which relates the 
difference between the spectral functions for vector and axial-vector currents to Am\ = m 2 K+ — 
m*p (DGMLY sum rule), which is a generalization of the Weinberg sum rule By assuming a 
saturation by the lowest resonance states of vector (rho meson) and axial-vector {a\ meson), their 
estimate of Am\ was already close to the experimental value, Ara|(Exp.) = 1261.2 MeV 2 [g], 
which implies that the non-perturbative contributions are important. There have also been results 
from extended chiral perturbation theory including resonance states [Q] or Bethe-Salpeter equation 
[Q]. In lattice QCD, Am\ has been calculated in quenched QCD with the Wilson fermion [^] and 
Nf = 2 QCD with domain- wall fermion []7|]. In these works, the EM interaction is introduced 
to make the (QCD+QED) system on the lattice, and the charged and neutral pseudoscalar meson 
masses are calculated through two-point functions in the usual way. Their values also show good 
agreement with Am|(Exp.) within the error. 

In this work, we apply the DGMLY sum rule to the evaluation of Am\ (for an early attempt, 
see Ref. With the DGMLY sum rule, A/n| is written in terms of a momentum integral of 

the difference of the vacuum polarizations defined by the vector and axial-vector currents. There 
are two points to be noted. First, the DGMLY sum rule exactly holds only in the chiral limit, and 
hence the lattice calculation requires good control of the chiral expansion. Second, the vector and 
axial-vector currents form a chiral multiplet in the continuum theory. Their difference signals the 
spontaneous breaking of chiral symmetry. With the domain-wall fermions the chiral symmetry is 
not good enough to calculate (VV — AA) unless the depth in the fifth dimension is unusually large. 
By using the overlap fermion, the chiral symmetry exactly holds, which makes the extraction of 
Am 2 n possible. 

The VV — AA vacuum polarization also provides pion decay constant f n and a low-energy 
constants (LECs) L r l(j . We also present a calculation of these quantities with two-flavor dynamical 
overlap fermion. 



2. Definition 



2.1 Continuum formula 

The leading order EM contribution to the pion mass difference is given by one photon ex- 
change diagram in the self-energy calculation of the charged and neutral pions as 

d 4 q 1 



Ami 



(lnf2 Dlxv{q) j d^liit+mj™,!™}^) - (*P\T{J™J™}\*P) 



(2.1) 



where D llv (q) is the photon propagator and = Y,f e fYfYnYf the EM current. Using the soft- 
pion relation and the current algebra, eq.([0|) can be written as 

" 2 ' d 4 q 



Ami 



f 2 J (2*) 4 



D, v (q) (o|r{^ 3 ,v v 3 }|o)(x)-(o|r{A3,A 3 v }|o)W 



(2.2) 
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3a 



EM 



*m Jo 



dQ 2 Q 2 n v _ A (Q z ), 



(2.3) 



where Q 2 = —q 2 >0,Vn = YYiiT a y and = ^ry^T"^ with SU(2) generator T a normalized 
by tiT"T b = 8 ab [jl|, g]. Hereafter f n = 130.5 MeV normalization is adopted. The second equation 
can be derived by substituting {0\T{J^J v }\0) = {d^Q 2 - QnQ v )Uj(Q 2 ) and D^ v = (8^ V Q 2 - 
(1 — ^)Q^iQv)/Q 4 with arbitrary gauge parameter E, into the first equation, and we define Ylv a = 
Ily — Ha- It should be noted that, in deriving eq. (jO|), the chiral limit is taken after applying the 
soft-pion relation. Thus the above is the exact formula in the chiral limit. 

CHPT at one-loop order predicts the low momentum behavior of Uv-a(Q 2 ) [10] as 



n v - A (Q 2 



f 2 

J 71 



Q 2 + ml 



8*.io(Mz) + 



1 



247T 2 L 



1 



-Iff (Jin - 

3 V a+1 



+ 2 



2 

-In-? 



+ 



(2-4) 

with a = \f\ + 4m 2 l /Q 2 . L r m is related to the S -parameter [11], which plays an important role in 
analyzing new physics models. On the other hand, at large momentum OPE [12] provides 



Cd= 2tn 2 (Q 2 ) C d=4 m q (yy) C^e 
Q 2 Q 4 Q 6 



+ 



(2.5) 



(2.6) 



(2.7) 



where the explicit form of Cd=2,4 has been known to two-loop order and Q/=6 is given by 

C d=6 = 87r(a,.<9 8 ) + (log term) + @{a 2 ), 

and 

(o 8 >* = ((n^°|vW W°~vO - m^ at ^w){wtys^w)) ^ 

A" is the Gell-Mann color matrix and t" is the Pauli matrix, and \i Q is a renormalization scale. The 
logarithmic term is estimated to be a few % of the leading term, which we ignore. (Og) contains 
information of the matrix element of ^° — > (nn)i = 2 [13]. 

2.2 Lattice formula 

In this simulation we use the following vector and axial-vector currents, 

A 



D m 



2mo 



D n 



2m Q 



(2.8) 



where mo = 1.6 and non-perturbative renormalization constant Zy = Za = 1.38 |[14|] is applied. 
Since these currents are not the conserved one, the current-current correlation functions may con- 
tain lattice artifacts. Their explicit form is represented by 



(0), 



+ (aQ ll fS^ v B 2 (Q i ) + (aQ^(aQ v y + (aQ^(aQ v ) )C n (Q z ) + 



(2.9) 



The ellipsis denotes the higher order terms and hereafter we ignore these terms. In order to re- 
move the lattice artifacts A,B\,B2,Cn we impose the Ward-Takahashi identity by constructing a 
set of linear equations, Q^{J^J V ) = fi(Tlj,A,B,C), Q^iQ v (J^iJ v ) = f2(Tlj,A,B,C), • ■ •, for differ- 
ent momentum configurations giving the same Q 2 . By solving these linear equations, we obtain 



n y = ni 0) +n 



i) 



3 



Pion mass difference from vacuum polarization 



E. Shintani 



-0.0001 



-0.0002 



-0.001 - 



-0.002 



-0.003 




0.04 0.06 0.08 

2 2 

a m 

71 



0.1 



Figure 1: The dependence of pion mass squared for Q 2 TIv-a at 4 low momenta, and 1 high momentum. 
The straight lines denote the fit function. Except for the lowest momentum fit function is linear function. 



3. Numerical results and analysis 



3.1 Lattice parameters 



We use the Nf = 2 dynamical overlap fermion configurations with the Iwasaki gauge action 
at /3 = 2.3 corresponding to aT l = 1.67 GeV on a 16 3 x 32 lattice generated at a fixed topolog- 



ical charge Q top = [15]. The sea quark masses are chosen to m q = 0.015, 0.025, 0.035, 0.050 
and the valence quark mass takes the same values as the sea. The pion mass squared are then 
0.082, 0.134, 0.189, 0.27 GeV 2 . The number of statistics is 200 configurations separated by 50 
HMC trajectories. The statistical error is estimated using the jackknife method with bin size equal 
to 2. 

3.2 Chiral extrapolation 

The chiral extrapolation of YIv-a(Q 2 ) is made at each momenta. The (a 2 m 2 ) dependence 
and fit results are shown in Fig.[T] for several representative (aQ) 2, s. At the lowest momentum, 



(aQ) = 0.0384, the CHPT prediction (2.4) is used together with f% and m 2 obtained from the 



pseudoscalar two-point correlation function. As shown in Fig.[j]the fit describes the data reasonably 
well, and leads to 

L r w (n = 770MeV) = -0.00474(23), (3.1) 

which is in good agreement with experimental value L' 10 (/i o = 770MeV)| Exp = —0.00509(47) 
[JT9]. For other momenta, the data show linear behavior, that we fit with a linear function. This im- 
plies that the quadratic quark mass term in eq.(|2~3|) is small, and the dimension-four and dimension- 
six terms are dominating in small quark mass region. 
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We also try to estimate the four-quark condensation (Og) by fitting to a functional form similar 

, _D 1 m 2 (Q 2 ) D 2 m q D 3 m.q + D 4 



to eq. Q2.5D, 



+ ■ 



■ + 



(3.2) 



Q2 Q4 Q6 

where the fit parameters D\, D2 and D^niq + D^ correspond to Cd=2, Q=4 (VV) an d Q=6> respec- 
tively, however any Q 2 dependence in Cd=2£fi is omitted. Fitting range is chosen as [1.235, 1.973] 
where OPE is expected to be a dominant contribution to the vacuum polarization. In the chiral 
limit, Da, = (8na s O&) up to a logarithmic correction, and we obtain 

-0.20(13) x 10~ 3 GeV 3 



2GeV 



(3.3) 



with a s (n = 2GeV) = 0.334, which is used in [Q]. In Ref.Q, their estimated value in the MS 
scheme is -(0.67 ~ 1.29) x 10~ 3 GeV 3 , which was the same sign with our result while whose 
magnitude is larger than eq.([0|). Although our result still contains a large uncertainty, it suggests 
a feasibility to extract the expectation values appealing in the QCD sum rule analysis. 

3.3 Numerical integral 

After chiral extrapolation we obtain the momentum dependence of Q 2 Yly-A in the chiral limit 
as shown in Fig.||. In order to perform the numerical integral for Q 2 H V A , we fit with an appropriate 
function. From eq. (2.4), Q 2 Yly A in the massless limit is given by 



lim Q 2 U V _ A (Q 2 ) = -f 2 n 



Q 2 . £P 

24k 2 h 2 



(3.4) 



which should be satisfied when making the fit ansatz. Furthermore, Yly a contains poles corre- 
sponding to resonance states at negative Q 2 . According to these requirements and the OPE predic- 
tion, we take the following fit function 



Q 2 nf_ A (Q z 



-F 2 + 



2 F 2 



Q 2 F : 



Q 2 ^Q 2 , n2 
-24k2~ +CiQ 



1 



(3.5) 



Q 2 +M 2 Q 2 +M 2 '\ 24n 2 ' J 1 +c 2 Q 6 ' 

A/12, F, F\ } 2,ci.2 are free parameters. The last term expresses the logarithmic term at small Q 2 and 
this is suppressed by C2Q 6 for large Q 2 . Figure || shows that the fit function (3.5) well describes 
our data below (ciQ) 2 = 2.0604. From the fitting result of F 2 , 



f„ = 107(15) MeV 



(3.6) 



is obtained. This value is consistent with f n 
hadron spectrum using the same configurations [1141]. 

In the numerical integral we split the integral range into two regions at Q 2 



Ami 



1 10 MeV, which is obtained from the study of 

A 2 : 

(3.7) 



-^TJlH dQ 2 Q 2 Uf_ A (Q 2 ) + f A2 dQ 2 Q 2 U^ A (Q 2 



Below A we use the fit function (3.5) and the parameters determined by the fit, while above A 



we employ the perturbative form at the one-loop order [ |12| ] with the factorization method for the 
expectation value of 0%, 



e 2 n° PE A 



64 n 



1 + 



K 



89 
48 



l, Q 2 
■In— ^ 



4~„2 



Q 



-4 



(3.8) 
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Figure 2: Momentum dependence of Q 2 TIv-a {Q 2 ) m the chiral limit. The dashed line denotes the prediction 
of OPE at leading order, and straight line denotes fit function. The dashed-dots line is a cutoff point. 

with ji = 2 GeV and (xj/yf) = -(251 MeV) 3 [[F7|]. A 2 is set to a 2 A 2 = 2.0604, where the gap 
between the two regions is negligibly small. Our result is 



Ami = 976(100) stat . +48 OPE MeV 2 = 1024(100) MeV 2 , 



(3.9) 



where we insert f n in eq.(|3.6[) into the denominator of eq. ( |3.7[ ). 

In addition to the statistical error given above, there are several sources of systematic error. 
Since the physical volume of our lattice is about (1.9 fm) 3 , the lightest pion data (~290 MeV) 
could receive sizable finite size effect. Also, since this result is obtained at a fixed topological 
sector (<2 top = 0), an additional finite size effect of 0(1 /V) is expected [19]. We may include such 



effects by modifying the chiral extrapolation. In the comparison to the experimental value, it is 
also necessary to evaluate the correction due to the small but finite quark masses. Although these 
systematic errors are yet to be estimated, it is encouraging that our result is reasonably consistent 
with the experimental value, Am 2 (Exp.) = 1261.2 MeV 2 . 



4. Summary 

In this work we have calculated the pion mass difference by applying the DGMLY sum rule 
to the lattice calculation for the first time. We also obtained f n and L r l0 by comparing the Q 2 de- 
pendence of the vacuum polarization with the predictions of CHPT. The use of the overlap fermion 
made this calculation possible. 
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